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We introduce a temperature-dependent parameterization in the modified embedded-atom method 
and combine it with molecular dynamics to simulate the diverse physical properties of the S and e 
phases of elemental plutonium. The aim of this temperature-dependent parameterization is to mimic 
the different magnitudes of correlation strength of the Pu 5/ electrons at different temperatures. 
Compared to previous temperature independent parameterization, our approach captures the neg- 
ative thermal expansion and temperature dependence of the bulk moduli in the 5-phase. We trace 
this improvement to a strong softening of phonons near the zone boundary and an increase of the /- 
like partial density and anharmonic effects induced by the temperature-dependent parameterization 
upon increasing temperature. Our study suggests it is important to include temperature-dependent 
parameterization in classical force field methods to simulate complex materials such as Pu. 



INTRODUCTION 

Plutonium, one of the most complex elemental met- 
als, has been a subject of intensive theoretical and ex- 
perimental research for more than five decades due to 
its unique anomalous physical properties (see [I] for a 
review). Plutonium has six equilibrium solid phases at 
ambient pressure but different temperatures. It under- 
goes a very large volume change (26 percent atomic vol- 
ume expansion) during the phase transition going di- 
rectly from the a phase to the S phase when the lat- 
ter is stabilized by a small concentration of impurities. 
The thermal expansion in the S phase of Pu is nega- 
tive, which is very unusual in metals; the melting point 
is very low relative to the other members of the ac- 
tinide series; the shear anisotropy of <5-Pu is the largest 
among all the elements; and the Gruneisen parameter 
of 5-plutonium is high compared with those of most ele- 
mental metals[2_. Furthermore, the elastic moduli of S- 
Pu have a surprisingly strong temperature dependence, 
softening with increasing temperature [5. The physical 
origin of the thermal properties, and in particular the 
negative thermal expansion in the S phase, remains mys- 
terious and has not yet been fully accounted for by first- 
principles electronic structure methods. The thermody- 
namics of elemental plutonium, while having some simi- 
larities with those of other heavy fermion materials, also 
pose unique challenges. The linear term of the electronic 
part of the specific heat is substantially enhanced over its 
local density approximation (LDA) band theory value. 
At higher temperatures, the temperature dependence of 
the electronic specific heat is anomalous. Different in- 
terpretations of this contributions as originating from 
the electron-phonon interaction and the electron-electron 
interaction have been discussed [IHS]. Furthermore the 
transition from a to i5 phase is characterized by a large 
excess of entropy beyond the phonon contribution, but no 
comprehensive microscopic theory of the excess entropy 
is available. [IllH] 

A very useful phenomenology for S-Pu and its alloys 
has been constructed by Lawson and his collaborators, 
by using analogies with the Invar model of iron |9l [10] . 



In this picture, Pu can be in two states which are as- 
sociated to very different volumes and energies, and the 
free energy can be used to reproduce the elastic and ther- 
modynamic properties of d Pu and their dependence on 
impurities. 

As pointed out by Johansson, [llj S-Pn sits at a spe- 
cial position in the actinide series, straddling between 
the so-called light actinides (Th to Np, with itinerant 
5/ electrons, smaller atomic volumes and metallic char- 
acteristics) and the heavy actinide elements (Am to Es, 
with larger atomic volumes and insulator-like localized 
5/ orbitals). The S-Pn f electrons are thus close to 
a localization-delocalization transition or crossover, in 
which their dual wave-like and particle-like character is 
important. ^12j This feature also appears clearly in the 
rearranged periodic table of Smith and Kmetko[T3]. 

Modern studies of this material, using a combination 
of density functional theory and dynamical mean field 
theory (DFT-hDMFT) have substantiated and refined 
this picture. [TiHTS] To describe plutonium, a large num- 
ber of atomic configurations are required, and valences 
with n/=5 and 6 have substantial weight in the ground 
state. This multi-determinantal mixed-valent character 
has been recently confirmed by resonant x-ray emission 
spectroscopy [12]. Unfortunately, anharmonic effects, and 
other effects which result from the motion of the atoms, 
are beyond the scope of current implementations of these 
first-principles techniques. They can be investigated us- 
ing powerful molecular dynamic (MD) methods, once a 
proper parameterization of the energy as a function of 
the atomic positions is posited. These methods require 
a parameterization of a classical many body potential. 
But with a small number of parameters, in a prescribed 
functional form, they can be used to predict successfully 
a large number of physical properties in a host of s and 
p band elements [^Dl - E5] . 

In an important series of papers, Baskes and 
collaborators |24f[26] used an extension of the highly- 
successful embedded-atom method (EAM) [171[1H], the 
modified as MEAM which includes explicitly angular 
terms and can describe more complex crystal struc- 
tures [20]. With this method, they described many 
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structural properties of elemental Pu at finite tempera- 
ture [24]. These calculations however do not reproduced 
the experimentally observed negative thermal expansion 
in ^-Pu. 

Force field methods such as MEAM can be in prin- 
ciple guided or parameterized using insights from first- 
principles calculations. The key part of a force field 
method is a formula E[Ri, R2, Rn] which expresses 
the total free energy in terms of fixed atomic positions 

of all the atoms i?i,i?2, ,Rn- It can be obtained by 

tracing out all the electronic coordinates, in a full path in- 
tegral treatment of ions and electrons, for a fixed position 
of the ions. In practice, for weakly correlated systems, 
it is approximated by a total energy calculations per- 
formed within the Born-Oppenheimer approximation at 
zero temperature but finite temperature corrections de- 
scribing the smearing of the Fermi function can be taken 
into account}^. 

In this work, we propose a simple modification 
of Baskes' parameterization of the MEAM force field 
method. We argue, that since in /-band materials like 
Pu, where strong Coulomb repulsions between electrons 
results in a strongly temperature dependent electronic 
structure, the many-body potential used in the MD sim- 
ulations has to contain an intrinsic temperature depen- 
dence. This temperature dependence in the form of 
Kondo resonance, is well established in cerium based ma- 
terials and is seen in the DMFT spectral function of many 
compounds, as for example elemental cerium and the 
heavy- fermion material Celrln^ [5D]. This temperature 
dependence has also been parameterized recently in the 
context of the two fiuid picture |31j . In Pu the character- 
istic energy scale T* is estimated to be about 800i4r[18]. 
Incorporating this idea into the MEAM parameterization 
of the free energy, we obtained with a small modification 
an improved description of the physical properties of Pu, 
including the negative thermal expansion of the S phase. 

MODEL AND COMPUTATIONAL DETAILS 

We use the atomistic model of plutonium formulated 
by Baskes |24j . except for introducing temperature de- 
pendence into the parameters /3*^°^ and (3^^\ At the heart 
of the MEAM formalism, the total energy of a system of 
monatomic atoms is given by 

Etot = J2ip(p^) + lJ2^i\^^M (1) 

i j^i 

where the first term on the right hand side is the em- 
bedding function F{pi), the energy required to embed an 
atom into the background host density pi at the lattice 
site i, which is calculated as 

F{p) = AE,^\n{^) (2) 

where A is an empirical parameter, po th^ background 
density of the reference structure and the cohesive 



energy. The second term, the pair potential 0, originat- 
ing from pairwise electrostatic repulsion between atoms 
at lattice sites i and j, is a function of the distance \rij\ 
between them. 

0(i?) = |[i?"(i?)-F(p(")(i?))] (3) 

where Z is the number of first-neighbor atoms in the 
structure and E^ is assumed the following form 

E'^iR) = -Se(l + a* +<5a*'^)e-''' (4) 
a* = 

with Tf. the nearest neighbor distance of the reference 
structure and fitting parameters a*, 5, and a. Back- 
ground electron density pi at site i is modeled as a combi- 
nation of angular contributions from neighboring atoms. 

p = p^o'yiTr (5) 

1=1 ^ 

where i*^'-' 2, 3) are fitting parameters, and the an- 

gular components of electron density p^'^ ('=0, 1, 2, 3) 
are calculated as weighted sum of atomic densities: 

= Y: P^aHRu) (6) 

where xf^ = Rf^/Rij, and pa^ are atomic electron densi- 
ties, taking an exponential form 

p(')(i?) = e-'''"W--i) (7) 

with fitting parameters /S*^'^ (^=0, 1, 2, 3). The second 
term in p^^^^ is added to make the partial electron densi- 
ties orthogonal. 

In the case of Pu, the atomistic model is parameterized 
by twelve semi-empirical MEAM parameters |24j. i.e., Ec, 
re, a, A, (1=0, 1, 2, 3), t^''> {1=1, 2, 3), and S, for 
which we take the same values as in Refill] except for 
two parameters discussed in details below. Among these 
parameters, of particular interests are (1=0-3), which 
characterize the exponential decay of s,p, d, /-like partial 
densities with the distance from the center of an atom. 
In RefI211 Baskes connected /3^'^ to the bonding char- 
acter of electronic orbitals. For Pu the most important 
parameters are Z?*^"-* which has the strongest effect on the 
potential energy and /S'-^-' which can mimic the changes of 
the 5/ electrons at different temperatures. As the emer- 
gence of all the phases of Pu is mainly due to thermal 



3 



effects on electronic correlations of 5/ electrons and the 
potential energies of the different phases, it is therefore 
natural to introduce a temperature-dependent parame- 
terization of P^^\T) and ^'^°^(T) to capture these effects. 

In this paper, we take /J^'^^ as a simple linear function of 
temperature, j3^^\T) = 9.3422 - This formulation 

closely follows the variation of the incoherent quasipar- 
ticle peak in 6-Pu, obtained from DMFT calculation [18]. 
For we also choose a simple linear function of tem- 
perature below 720 K and fix it for temperature above 
720 K, considering the fact that Pu undergoes the S-e 
phase transition around 720 K which temperature is also 
close to the coherence temperature of Pu[18]. Specifi- 
cally we choose I3(-"^T) = 1.634 + 0.0014 * T for T<720 
K and l3^°^iT) = 2.642 for T>720 K. We show that such 
a temperature dependent parameterization can account 
for several experimental observations which cannot be 
described by the previous temperature independent pa- 
rameterization. 

Compared to the previous temperature independent 
study[23], we include many more atoms and allow much 
more time for the system to achieve better thermal- 
dynamical equilibrium. Specifically, we use 6x6x6 su- 
percell (864 atoms) for the 5 phase and 8x8x8 supercell 
(1024 atoms) for the e phase. We run the MD simulation 
for 450 ps to allow the system to arrive at good thermal- 
dynamical equilibrium and use the next 50 ps to compute 
the average values of relevant physical quantities such as 
volume. 

RESULTS AND DISCUSSION 

Using the proposed temperature-dependent parame- 
terization, we performed the MEAM-I-MD calculations 
to compute several physical properties, including the vol- 
ume, bulk moduli, and phonon density of states. 

Volume 

We first show in Fig{l] the calculated volumes of S and 
e plutonium at different temperatures using the tempera- 
ture dependent parameterization and compare them with 
the corresponding results of the previous calculations ^i] 
using temperature independent parameterization. Also 
shown in Fig|l] are the experimental values taken from 
Ref. 1321 As mentioned before, an important hallmark 
of the S phase is its negative thermal expansion, i.e., its 
volume decreases with increasing temperature. This un- 
usual phenomenon doesn't yet receive a satisfactory the- 
oretical explanation. Clearly, the calculated volume us- 
ing the temperature independent parameterization has 
a very weak temperature dependence in the S phase, 
and this negative thermal expansion is not captured by 
the temperature independent parameterization. However 
with our temperature dependent parameterization, we re- 
cover the negative thermal expansion in the S phase. This 
improvement suggests that thermal effects on the 5/- 
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FIG. 1: (color online) The calculated volumes of S and e 
plutonium at different temperatures using the temperature 
dependent parameterization (this work) and temperature in- 
dependent parameterization (data taken from Ref. I24p and 
compared with experimental values taken from Ref. 1321 

electronic correlations must be effectively incorporated 
into the force field method in order to recover the phe- 
nomenon of negative thermal expansion. For the e phase, 
the calculated volume by the temperature independent 
parameterization is almost fiat as a function of temper- 
ature, whereas our temperature dependent parameteri- 
zation captures the positive thermal expansion in this 
phase, thus improves again previous results of the tem- 
perature independent parameterization. 

Bulk modulus 
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FIG. 2: (color online) Comparison of the calculated Bulk 
modulus in 5-Pn using temperature-independent and our pro- 
posed temperature-dependent parameterizations and with ex- 
perimental data taken from RefjS] 
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We now proceed to compare the bulk moduli from the 
temperature dependent and independent parameteriza- 
tions to experimental studies. In recent years, resonant- 
ultrasound spectroscopy measurements have been able to 
study the elastic moduli of various phases of Pu over a 
wide range of temperatures [21 131 133 [311 ■ The most re- 
markable common observation from these experiments is 
the large elastic softening of the bulk and shear mod- 
uli with increasing temperature. Such unexpected an- 
harmonic behavior in plutonium has been conjectured 
to arise due to 5/-electrons undergoing localization- 
delocalization transition. In Figure. [2| we compare the 
calculated bulk moduli using our proposed temperature- 
dependent parameterization with that using traditional 
temperature-independent one for the S phase. Experi- 
mentally, the bulk modulus was measured up to about 
SOOi^T, in the stabilized fee structure [3], which shows 
monotonic decreasing B{T) as a function of temperature. 
Although the measurement was not extended to higher 
temperature, the softening tread is expected to continue. 
In our MD calculations, the temperature-dependent pa- 
rameterization produces the elastic softening with in- 
creasing temperature, consistent with experiments. This 
comparison shows that the observed softening of elastic 
moduli originates from the 5/ electronic correlations in 
Pu, rather than a lattice effect. 

Phonon density of states 

We also explore the phonon density of states (DOS) in 
5 plutonium. Experimentally, inelastic neutron scatter- 
ing (INS) experiments by McQueeney et al. [3 [3S] stud- 
ied the temperature dependence of the phonon spectrum 
as observed from measurements of phonon density of 
states (DOS). The studies revealed unusual lattice behav- 
ior, namely a decrease in phonon energies with increas- 
ing temperature, which again was ascribed to effects of 
strong electronic correlations of the Pu 5/ electrons. The 
phonon DOS was not presented in the early MEAM-t-MD 
study of Pu [21] nor was shown in a following lattice vi- 
bration study of the S-Pu [25] using the MEAM-hMD 
method. Here in the top panel of Fig. [3j we present 
calculated phonon density of states at two different tem- 
peratures {580K and 720K ) both with and without the 
temperature-dependent parameterization. In the bottom 
panel of Fig. [3] we compare our temperature-dependent 
calculation with experimental findings [71 135] and previ- 
ous theoretical work of first-principles calculations us- 
ing DFT-HDMFT[Tn]. Ahhough it is difficult to see the 
phonon softening effect in our calculated phonon DOS re- 
sults, the main features of our calculated density of states 
is in good agreement with the measured phonon density 
of states. This agreement verifies the correct physics rep- 
resented by the MEAM parameters, including our new 
temperature-dependent ones. The MD calculations cap- 
ture the low frequency peak at around d.^meV, which 
is not visible in the first-principles calculation T5]. How- 
ever our calculated phonon DOS are slightly off at higher 
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FIG. 3: (color online) Top panel: Calculated phonon den- 
sity of states in S-Pu with temperature independent and 
temperature dependent parameters. Bottom panel: Com- 
parison of phonon density of states in 5-Pu with experi- 
mental measurements [7] and first-principles theoretical 
calculations [iSj . 

frequencies. This discrepancy may be due to tempera- 
ture effect since our calculations in the 6-Pu are at 580K 
whereas experimental measurements are at room temper- 
ature, and in addition the experiments have some impuri- 
ties in the delta-Pu (5 percent Al in Ref.[35land 2 percent 
Ga in Rcf. 7) In the following section, we show that using 
the MEAM formalism alone without the MD simulations, 
our temperature-dependent parameterization but not the 
temperature independent parameterization captures the 
phonon softening effects in the S phase. 

Qualitative understanding 

To understand the physical consequences introduced 
by the temperature-dependent MEAM parameters (/?'-*'■* 
and Z?*-"^^ in this work), we have calculated the phonon 
dispersion of S-Pu within the harmonic approximation, 
where the dynamical matrices are evaluated solely from 
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FIG. 4: Phonon dispersion and phonon density of states of 
S-Pu with lattice constant a — 4.64^4, calculated analytically 
from the MEAM model within the harmonic approximation. 
The temperature dependence is incorporated by using the 
temperature dependent /3'°>(r) and /3(^'(r). The phonons 
near the zone boundary soften strongly upon increasing tem- 
perature, which is consistent with the negative thermal ex- 
pansion. 



the total energy given by the MEAM formalism without 
an explicit temperature dependence as included in the 
MD simulations. Here the temperature dependence is in- 
corporated by using the temperature dependent /3'^*'^(T) 
and /3(3)(T). 

Our results are shown in Fig[4j An surprising observa- 
tion is that upon increasing temperature, clear and dras- 
tic phonon softenings occur at the X and L points indi- 



cating that the i5-phase tends to be unstable with increas- 
ing temperature, in reminiscent of the unstable phonons 
of e-Pu in DMFT calculations [E] by Dai et al. Mean- 
while, the overall scale of the phonon spectra, phonon 
density of states and the elastic constants (Cii,Ci2,C4, 
derived from the long-wave phonon dispersion) change 
only slightly, consistent with the results given by the MD 
simulations. The long-wave phonon dispersion agrees 
quite well with the experimental observation that the TA 
and LA branches are nearly degenerated along the T-X 
direction. Also the elastic constants extracted from 
the dispersion agree well with experimental results. [57] 
These agreements justify the validity of the current tem- 
perature dependent parameterization. On the other 
hand, within Baskes's temperature-independent param- 
eters, the phonon dispersion has no temperature depen- 
dence within the harmonic approximation (the only pos- 
sible temperature dependence is from the change of the 
lattice constants against the change of temperature, how- 
ever the lattice constant of the temperature independent 
parameterization has a very weak temperature depen- 
dence as shown in FigjlJ. Thus, the main effect intro- 
duced by our temperature-dependent parameter is the 
inclusion of phonon-softening effects at the level of the 
harmonic approximation. These effects seem to be es- 
sential to obtain the correct negative thermal expansion 
in the (5-phase. 

The results also indicate that anharmonic effects in- 
duced by thermal vibrations (molecular dynamics) play 
an important role in stabilizing the 5 phase in the MD 
simulations. The phonon softening in the harmonic ap- 
proximation is so drastical that the J-phase should be 
unstable when temperature is higher than 580 K where 
the (5-phase is actually stabilized. This implies that the 
(5-phase is stabilized by the MD procedure with the help 
of anharmonic effects in phonons. This could also be the 
possible reason why at the zone boundary the calculated 
phonon dispersions do not agree well with experimental 
results and also why the overall scale of phonon spectra 
(around 10 meV in the calculation) is about 20% smaller 
than in the experiments [71 [55j and the MD simulations. 

For cubic crystal structures such as the 5 and e plu- 
tonium, the anharmonic effects introduced in the MD 
simulations are reflected in the non-s angular partial den- 
sities /O*^'-* (/ — p, d, /) which are zero in the MEAM meth- 
ods and acquire finite values during the MD simulations. 
These non-s angular partial densities p^'^ (/ = p, d, /) 
characterizes the the weight of the non-spherically sym- 
metric fluctuations in the positions of the neighboring 
atoms. The bigger the non-s angular partial densities, 
the larger the anharmonic effects introduced and the de- 
viations from the cubic crystal structures. Therefore, the 
non-s angular partial densities can also be viewed as a 
measurement of the randomness of the atoms away from 
the cubic structures during the MD simulations. We find 
that, compared to the temperature-independent parame- 
terization, the temperature-dependent parameterization 
has more /-angular-dependence in the partial charge den- 
sity in 5 Pu. This suggests that it is important to include 
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appropriate amount of /-angular-dependent partial den- 
sity in 6 Pu in order to account for experimental obser- 
vations, and introducing the temperature-dependent /?'^-' 
and is an efficient way to fulfill this purpose. 

While invoking our proposed temperature-dependent 
parameterization, we argue that the conventional 
temperature-independent parameterization misses some 
important contributions from the 5/-electrons in the high 
temperature (e.g. 6) phases, and hence it fails to describe 
certain physical properties (like negative thermal expan- 
sion) in the (5— phase, which arise due to strong corre- 
lations of 5/ electrons. Our temperature-dependent pa- 
rameterization introduces strong phonon softening upon 
increasing temperature and enhances anharmonic effects 
which help to stabilize the structure in the MD simu- 
lations. A combination of both effects accounts for the 
unique features of the negative thermal expansion and 
the phonon softenings in the S-Pu. 



CONCLUSION 

In this paper, we have demonstrated that an intrin- 
sic temperature dependent force field method is needed 
to simulate the different phases of Pu, especially the 
S phase. To incorporate this temperature dependence, 
we use temperature dependent /3'^' and /S^'^-' parameters 
in the MEAM-I-MD calculations. The negative thermal 
expansion and the temperature dependence of the bulk 
modulus in the (5-phase, as well as the positive thermal 
expansion in the e phase are improved over the original 
temperature- independent MEAM-I-MD results and are in 
reasonable agreements with experiments. Our temper- 
ature dependence parameterization provides the simple 
idea that the large changes in the spectral function of Pu 
that were shown to take place in the electronic structure 
calculations [18] should affect the effective force between 
the Pu atoms which are described by the MEAM force 
field. The increase in upon increasing temperature 
incorporates several important effects compared to the 
temperature independent parameterization, (1) it intro- 
duces strong phonon softenings enhancing anharmonic- 
ity; (2) it enhances the weight of the non-spherically 
symmetric fluctuations in the positions of the neighbor- 
ing atoms; (3) it lowers the energy of the competing 
phases which have smaller volumes. Fluctuations into 
those phases decrease the volume in a way reminiscent of 
the Invar model[9l[T^, without having to include explicit 
Ising degrees of freedom as was done recently by Lee 
et al. in Ref.l^ Our results demonstrate the require- 
ment for finding more accurate atomistic potentials with 
temperature-dependent parameterization to incorporate 
the thermal effects on the electronic contributions to the 
system. To put the parameterization on solid theoretical 
grounds, there is the need for a microscopic derivation 
of the atomistic potentials including the temperature de- 
pendence starting from first-principles electronic struc- 
ture methods. Recent advances in the implementations 
of DFT+DMFT methodologies (see RefsOHlllDlfor a cou- 



ple of examples) look very promising in this respect. 



ACKNOWLEDGEMENTS 

This work was supported by DOE BES grant num- 
ber DE-FG02-99ER45761. K.B. wishes to thank the 
Alcatel-Lucent foundation for support. Q.Y. was sup- 
ported by the DOE nuclear Energy University Program, 
contract No. 00088708. The authors generously thank 
Mike Baskes and Steve Valone for sharing their MD code 
and for helpful discussions. 



[8 

[9 
[10 

[11 
[12 

[13 
[14 
[15- 
[16: 
[17 
[18 
[19 



[20" 
[21 



N. G. V. Ed., Challenges in Plutonium Science, Los 
Alamos Science 26 (2000). 

H. Ledbetter, A. Migliori, J. Betts, S. Harrington, and 
S. El-Khatib, Phys. Rev. B 71, 172101 (2005). 
A. Migliori, H. Ledbetter, A. C. Lawson, A. P. Ramirez, 
D. A. Miller, J. B. Betts, M. Ramos, and J. C. Lashley, 
Phys. Rev. B 73, 052101 (2006). 

M. J. Graf, T. Lookman, J. M. Wills, D. C. Wallace, and 
J. C. Lashley, Phys. Rev. B 72, 045135 (2005). 
J. C. Lashley, J. Singleton, A. Migliori, J. B. Betts, R. 
A. Fisher, J. L. Smith, and R. J. McQueeney, Phys. Rev. 
Lett. 91, 205901 (2003). 

A. C. Lawson and J. C. Lashley, to be published. 
M. E. Manley, A. H. Said, M. J. Fluss, M. Wall, J. 
C. Lashley, A. Alatas, K. T. Moore, and Yu. Shvyd'ko, 
Phys. Rev. B 79, 052301 (2009). 

J. R. Jeffries, M. E. Manley, M. A. Wall, K. J. M. 
Blobaum, and A. J. Schwartz, Phys. Rev. B 85, 224104 
(2012). 

A. C. Lawson, J. A. Roberts, B. Martinez, and J. W. 
Richardson, Philos. Mag. B 82 1837-1845 (2002). 

A. C. Lawson, J. A. Roberts, B. Martinez, M. Ramos, G. 
Kothar, F. W. Trouw, M. R. Fitzsimmons, M. P. Hehlen, 
J. C. Lashley, H. Ledbetter, R. J. Mcqueeney, and A. 
Migliori, Philos. Mag. 86 2713-2733 (2006). 

B. Johansson, J. Magn. Magn. Mater. 47-48, 231 (1985). 
G. H. Lander, Science 301, 1057-1059 (2003). 
J. L. Smith and E. A. Kmetko, J. Less Common Met 90, 
83-88 (1983). 

S. Y. Savrasov, G. Kotliar, and E. Abrahams, Nature 
410, 793-795 (2001). 

X. Dai, S. Y. Savrasov, G. Kothar, A. Migliori, H. Led- 
better, and E. Abrahams, Science 300, 953-955 (2003). 

A. Toropova, C. A. Marianetti, K. Haule, and G. Kotliar, 
Phys. Rev. B 76, 155126 (2007). 

J. H. Shim, K. Haule, and G. Kothar, Nature 446, 513- 
516 (2007). 

C. A. Marianetti, K. Haule, G. Kotliar, and M. J. Fluss, 
arXiv:0805.1604vl. 

C. H. Booth, Y. Jiang, D.L. Wang, J.N. Mitchell, PH. 
Tobash, E.D. Bauer, M.A. Wall, P.G. Allen, D. Sokaras, 

D. Nordlund, T.-C. Weng, M.A. Torrez, and J.L. Sarrao 
PNAS 109, 10205-10209 (2012). 
M. I. Baskes, Phys. Rev. B 46, 2727-2742 (1992). 

B. J. Lee, J. H. Shim, and M. I. Baskes, Phys. Rev. B 
68, 144112 (2003). 

[22] Y.-M. Kim, B.-J. Lee, and M. L Baskes, Phys. Rev. B 
74, 014101 (2006). 



7 



[23] B. Jclinck, J. Houze, S. Kim, M. F. Horstcmcycr, M. I. 

Baskcs, and S.-G. Kim, Phys. Rev. B 75, 054106 (2007). 
[24] M. I. Baskcs, Phys. Rev. B 62, 15532-15537 (2000). 
[25] M. I. Baskes, A. C. Lawson, and S. M. Valone , Phys. 

Rev. B 72, 014129 (2005). 
[26] S.M. Valone, M.I. Baskcs, and S.P. Rudin, Journal of 

Nuclear Materials 422, 20 (2012). 
[27] M. S. Daw and M. I. Baskes, Phys. Rev. Lett. 50, 1285- 

1288 (1983). 

[28] M. S. Daw and M. I. Baskes, Phys. Rev. B 29, 6443-6453 
(1984). 

[29] G.J. Ackland, arXiv:1205.6684. 

[30] J. H. Shim, K. Haule, and G. Kotliar, Science 318, 1615- 

1617 (2007). 

[31] Y.-F. Yang and D. Pines, Phys. Rev. Lett. 100, 096404 
(2008). 

[32] O. J. Wick, Plutonium Handbook: A Guide to the Tech- 
nology, Vol. I and II (The American Nuclear Society, 
1980). 

[33] A. Migliori, C. Pantea, H. Lcdbcttcr, I. Stroc, J. B. Bctts, 
J. N. Mitchell, M. Ramos, F. Freibert, D. Dooley, S. 



Harrington, and C. H. Mielke, J. Acoust. Soc. Am. 122, 

1994-2001 (2007). 
[34] L. Bourgeois, M.-H. Nadal, F. Clement, G. Ravel- 

Chapuis, J. Alloys Compd. 444-445, 261-264 (2007). 
[35] R. J. McQueeney, A. C. Lawson, A. Mighori, T. M. Kel- 

ley, B. Fultz, M. Ramos, B. Martinez, J. C. Lashley, and 

S. C. Vogel, Phys. Rev. Lett. 92, 146401 (2004). 
[36] J. Wong, M. Krisch, D. L. Farbor, F. OcceUi, A. J. 

Schwartz, T.-C. Chiang, M. Wall, C. Boro, and R. Xu, 

Science 301, 1078-1080 (2003). 
[37] H. M. Ledbetter and R. L. Moment, Acta Metallurgica 

24, 891 (1976). 

[38] T. Lcc, M. I. Baskes, A. C. Lawson, S. P. Chen, and S. 

M. Valone, Materials 5, 1040-1054 (2012). 
[39] K. Haulc, C.-H. Yee, and K. Kim, Phys. Rev. B 81, 

195107 (2010). 

[40] M. Aichhorn, L. Pourovskii, V. Vildosola, M. Ferrero, 
O. ParcoUet, T. Miyake, A. Georges, and S. Biermann, 
Phys. Rev. B 80, 085101 (2009). 



